Description of the RHIC p^-spectra in a thermal model with expansion! 
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The assumption of simultaneous chemical and thermal freeze-outs of the hadron gas leads to a 
surprisingly accurate, albeit entirely conventional, explanation of the recently-measured RHIC p±- 
spectra. The original thermal spectra are supplied with secondaries from cascade decays of all 
resonances, and subsequently folded with a suitably parameterized expansion involving longitudinal 
and transverse flow. The predictions of this thermal approach, with various parametrizations for 
the expansion, are in a striking quantitative agreement with the data in the whole available range 
of < p_L < 3.5GeV. 
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In this Letter we offer a very simple explanation of 
the ^)_L-spectra recently measured at RHIC Our 
approach has the following ingredients: i) simultaneous 
chemical and thermal freeze-outs, with the hadron dis- 
tributions given by the thermal model; in other words, 
hadrons decouple completely when the thermodynamic 
parameters reach the freezing conditions, and no particle 
rescattering after freeze-out is present, ii) these thermal 
distributions are folded with a suitably parameterized hy- 
drodynamic expansion, involving longitudinal and trans- 
verse flow, finally, iii) feeding from resonances, including 
cascades, is incorporated in a complete way. 

So far, the thermal approach has been applied suc- 
cessfully in studies of particle ratios measured in rel- 
ativistic heavy-ion collisions at AGS and SPS [p|-p^. 
Quite recently, it has also been shown that particle ra- 
tios measured at RHIC may be equally well described 
in the framework of such models [pT|-p^ . Description of 
hadronic p_L-spectra in thermal models is more involved, 
since the spectra are affected by decays of resonances, 
hydrodynamic flow, and possibly by other phenomena 
occurring during the alleged phase transition from the 
quark-gluon plasma to a hadron gas |l4j. The model 
results presented in this Letter are in a surprising agree- 
ment with the experiment in the entire range of the data, 
< _p± < 3.5GeV, as can be seen in Fig. 1. The model 
has two free parameters: one controlling the size of the 
system (overall normalization of the spectra), and the 
other one the transverse flow. We test two different mod- 
els (parametrizations) for the freeze-out hypersurface and 
the hydrodynamic expansion. Both combine the Bjorken 
expansion with transverse flow , and follow the 

spirit of Refs. ||-||. 

The first model (model I) assumes that the freeze- 
out takes place at a fixed value of the invariant time, 

T = ^t'^ ~ ~ ''k ~ = const, which means that, due 

to time dilation, the particles in the fluid elements mov- 
ing farther away from the collision center decouple later 
than the particles in the fluid elements remaining at rest 
in the center-of-mass system of the colliding nuclei. Fur- 
thermore, we assume that the four- velocity of expansion 



is proportional to the coordinate, 
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The freeze-out hypersurface is parameterized as 
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where ay is the rapidity of the fluid element (vz — r^jt = 
tanhQ!||), whereas describes the transverse size of the 
system (p = rsha^). The transverse velocity is Vp = 
tanh a^/cha|| . We account for the finite transverse size 

of the system by imposing the condition p - 



n < 



Pmax- The second model considered (model II) has the 
same parametrization of the four- velocity and the freeze- 
out hypersurface as the popular blast model [19 2l|]: 



\^ = (cha||ch/3x, sha||ch/3j_, cos0sh/3x, sin0sh/3x) , 
t ~ rchaii , 
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Both models I and II are boost-invariant. Since devia- 
tions from boost-invariance are seen in the rapidity distri- 
butions at RHIC 1^, our present approach should be re- 
garded as an approximate treatment of the mid-rapidity 
region. However, this approximation is very good. One 
may depart from the boost-invariance by limiting the in- 
tegration in the a\\ variable, thus limiting the longitudi- 
nal size of the system. The results for the spectra ob- 
tained that way are very similar to the boost-invariant 
results presented below even for the limit for a\\ as low 
as 0.5. 

The local freeze-out conditions, i.e., the values of the 
temperature and the chemical potentials, are universal 
for the whole freeze-out hypersurface. Since for boost- 
invariant models the particle ratios at mid-rapidity are 
not affected by the expansion (this important point is 
discussed below), the values of the thermodynamic pa- 
rameters may be obtained directly from the standard 
thermal analysis, which yields |l3| T = 165 ± TMeV, 
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fiB = 41 ± 5MeV, fis = 9MeV, and ^/ = -IMeV. 
Knowing T and /i's we calculate the local distribution 
functions of hadrons which include the initial thermal 
contribution, as well as additional contributions from the 
sequential two- and three-body decays of all heavier res- 
onances. These decays are very important, since they 
effectively cool the system by 35-40 MeV, as recently 
shown in Ref. p3[ |, and also known from earlier works 
on other reactions p3| , p^ . The standard Cooper- Frye- 
Schonberg formula ]25[ | is used to calculate the p_L-spectra 
of the observed hadrons. As the result, the particle 
densities are obtained as the integrals over the freeze- 
out hypersurface, Ni = J d^p/p'^ J p^d'E^fi{p ■ u), where 
dT,fj_ is the volume element of the hypersurface, fi is 
the phase-space distribution function for particle species 
i (composed from the initial and secondary particles), 
and p^ = {m ±chy, m ±shy, p± cos Lp,p± sin ip) is the four- 
momentum. Finally, we find, for the case of model I, 
the rapidity and transverse-momentum distributions of 
hadrons, 
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where p ■ u — m^ch (j; — q:[|) cha^ — cos^shaj_ and 
^ — (f> ~ (p. Similarly, for model II the relevant parame- 
ters are the velocity, P±, and the volume of the system 
1^,0. We observe that the rapidity distribution (^) 
is boost-invariant, since the dependence on y can be ab- 
sorbed in the integration variable by shifting an — > Q!|| -\-y. 
Clearly, this is a direct consequence of the assumed boost- 
invariant form of the freeze-out surface. 

In the thermal-model fit of particle ratios o f Ref. fl^ ] 
one computes the integrals Ni = V J d^p fi{\/mf 
The question arises whether the ratios obtained that way, 
which correspond to the collection of particles from the 
whole phase space, are the same as the ratios experimen- 
tally measured in the mid-rapidity region, i.e., the ratios 
of the integrals dNi/dy = / d'^p±dNi/{d'^p±dy). The an- 
swer is yes, and follows from the boost-invariance of the 
expansion model. Indeed, since dNi/dy are independent 
of 2/, we have 



dNj/dy ^ J dydNjdy _ Ni_ 
dNj/dy ~ J dy dNj/dy ~ Nj 
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This obvious general result can be verified explicitly in 
our specific boost-invariant models. 

Figure 1 shows our main result. In Fig. 1 (a,b) we 
compare the model predictions for the p^-spectra of pi- 
ous, kaons, protons and antiprotons, with the PHENIX 
minimum bias preliminary data 0]. The model param- 
eters are fitted by the method including all points 
in the range < p_L < 2GeV from Fig. 1 (a,b). This 



yields r = 5.55fm, pmax = 4.50fm for model I, and 
V = (6.48fm)3, /3_L = 0.52 for model II. The statisti- 
cal errors in the fitted parameters are of the order of 
1%. In Fig. 1 (c) we compare the model predictions 
to preliminary PHENIX |§ and STAR ^ data for the 
most central collisions. The model parameters are fit- 
ted by the least-square method to all points in Fig. 1 
(c), yielding r = 7.66fm, pmax = 6.69fm for model I, 
and V = (9.81fm)3, = 0.52 for model II. It is natu- 
ral that the size parameters of the system for the most 
central collisions are larger than for the minimum-bias 
case, which averages over centralities. The size param- 
eters correspond to the size of the system, in particular 
for central collisions in model I p,nax = 6.69fm is very 
close to the radius of the Au nucleus, 6.22fm. We note 
that the quality of the fit in Fig. 1 is impressive. For the 
minimum-bias data (Fig. 1 (a,b)) model I (thicker lines) 
crosses virtually all data points within error bars, with 

/degree of freedom less than 1. Amusingly, also the 
very-high p± data are reproduced. The fit with model 
II is very similar in the range < p_L < 2GeV, and falls 
slightly below the data at higher p±, where hard pro- 
cesses are expected to contribute. The fit to the most 
central collisions. Fig 1 (c), is of similar quality except 
for the p preliminary data from STAR, which fall 30-50% 
below the model fits at low p±. Note, however, that at 
low p± the STAR and PHENIX preliminary data are not 
fully consistent. The total multiplicity of charged par- 
ticles produced in the most central event measured at 
PHOBOS [|7), 555 ± 12 ± 35, is consistent with size pa- 
rameters quoted above, and yields r 2± 7fm for model I 
and V ~ (9fm)'^ for model II. 

Since the values of the strange and isospin chemical 
potentials are very close to zero, the model predictions 
for 7r+ and 7r~, as well as for K'^ and K~ are practically 
the same. The value of the baryon chemical potential of 
41 MeV splits the p and p spectra. Note the convex shape 
of the pion spectra in Fig. 1, reproduced by the model. 
In addition, the 7r+ and p curves in Fig. 1 cross at p± ~ 2 
GeV, and the A'+ and p at ~ 1 GeV, exactly as in 
the experiment. We stress that our method is different 
from traditional fits in the blast or similar models, where 
the temperatures for various particles are being adjusted 
independently or from semi-empirical formulas. We have 
no freedom here: the temperature is the freeze-out tem- 
perature fixed by the particle ratios, and the spectra are 
obtained as described in Ref. [|3j. 

In Table 1 we present the inverse slope parameters, 
Tcff, defined by fitting the function const exp(—mj^/Tcff) 
to the data |^, and average p±. Agreement of model I 
and the data for the inverse slopes is within error bars 



* Note that this definition depends on the fitting region in 
m±. Thus, the resulting numbers depend quite strongly on 
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except for most central p data from PHENIX (2.5 stan- 
dard deviations) and K~ and p data from STAR (2 and 
4.6 standard deviations, respectively). The values of {p±) 
are within error bars. 

We end this Letter with a more pedagogical discus- 
sion of the role of various effects included in our anal- 
ysis. In Fig. 2 the dotted line shows the initial pion 
p_L-spectrum in a static fireball with the same temper- 
ature and chemical potentials as used in our calcula- 
tion. No secondaries are included here. The effect of 
the decays of all resonances is represented by the dashed- 
dotted line, representing the sum of the initial and sec- 
ondary pions. Decays of the resonances lead to an ef- 
fective decrease of the temperature by about 35-40 MeV 
1^ , since the emitted particles tend to populate the low- 
p± region. However, the spectrum remains concave. The 
effect of the pure longitudinal Bjorken expansion (with 
T = ^t"^ — = const) is illustrated by the dashed line. 
This is a redshift effect, since all fluid elements move away 
from the observer, which leads to extra cooling of the 
spectrum. The solid line corresponds to model I (model 
II gives very similar results), incorporating both the lon- 
gitudinal expansion and the transverse flow. The trans- 
verse flow causes some fluid element to move in the di- 
rection of the observer, leading to blueshift Hence 
we find a combination of redshift and blueshift, yield- 
ing the p± spectrum displayed by the solid line in Fig. 
3. Note that the spectrum finally acquires the convex 
shape, as seen in the experiment (c/. Fig. 1). The effects 
of blueshift are stronger for more massive particles, hence 
the behavior of Fig. 1 and Table 1. 

To conclude, we emphasize that the presented descrip- 
tion, implementing in a simple fashion all key ingredients: 
freeze-out, decays of resonances, and longitudinal and 
transverse flows, works for RHIC in the whole available 
range of data. We stress that we have assumed that the 
chemical and thermal freeze-outs occur simultaneously, 
which is in the spirit of the sudden hadronization of Refs. 
p9| , ^ . A practical value of our results is that they give 
hints and constraints for more involved hydrodynamic 
calculations, e.g. [^-^, by providing the freeze-out 
conditions that describe the data. A natural extensions 
of the model should include different centrality effects 
(elliptic flow), and rapidity dependence. Also, the model 
must be further verified against the available data from 
the HBT pion interferometry |Q . Note that our size pa- 
rameters are very similar to the experimental HBT radii. 
We have checked that our model works also for the SPS 



this region, and Tch is a largely-biased measure of the particle 
spectra. Following Ref. [§, we use O.SGeV < p± < 0.9GeV 
for the pions, 0.55GeV < p± < 1.6GeV for K'^ , p, and p, and 
0.75GeV <p± < 1.6GeV for K~ . The fits are performed 
in these regions. 



data. The details of these studies will be presented else- 
where. 
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FIG. 1. The px-spectra at j/ = of pions (solid line), kaons (dashed line) and protons or antiprotons (dashed-dotted line), 
as evaluated from model I (thicker lines) and model II (thinner lines), compared to the PHENIX preliminary minimum-bias 
data (a,b), and to STAR (open symbols) and PHENIX preliminary highest-centrality data (c) {Au + Au at 130GeV ). 
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FIG. 2. Contributions of various effects to the p^-spectra 
of Tv^ (normalizations arbitrary). 



TABLE I. Inverse slope parameters, T^b, and the average 
transverse momentum, (p±). 



Toff [MeV] 


TT 


K 


p or p 


PHENIX, min. bias, exp. 


190 ± 17 


272 ± 16 


274 ± 11 


positive hadrons model I 


203 


263 


292 


PHENIX, min. bias, exp. 


207 ± 20 


260 ± 17 


301 ± 14 


negative hadrons model I 


202 


255 


291 


PHENIX, most central, exp. 


197 ± 19 




367 ± 23 


negative hadrons model I 


206 


273 


310 


STAR, most central, exp. 


188 ± 20 


300 ± 35 


560 ± 50 


negative hadrons model I 


176 


223 


327 


(px> [MeV] 


TT 


K 


p or p 


PHENIX, most central, exp. 


370 ± 70 


600 ± 60 


840 ± 50 


positive hadrons model I 


438 


629 


853 


PHENIX, most central, exp. 


370 ± 70 


630 ± 80 


860 ± 50 


negative hadrons model I 


434 


629 


852 
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